ASYMPTOTIC PROPERTIES OF THE COLUMNS IN THE 
PRODUCTS OF NONNEGATIVE MATRICES 



ERIC OLIVIER AND ALAIN THOMAS 



Abstract. We consider the sequence of column vectors Ai . . . AnV associated to a se- 
quence (A„)„gN of nonnegative d x d matrices and to a positive d-dimensional column 
vector V. We give some sufficient conditions on (A„)„gN for the sequence of normal- 
ized column vectors -r—-^ — - to converge (do not confuse with the convergence of 

\\Ai...AnV\\ 

Al An 

-— — that can fail even in case there exists two positive matrices Mi and M2 

\\Ai . ..A„\\ 

such that An G {Mi,M2} for any n). These conditions seem to be strong but, given a 
sequence (A„)„gN of nonnegative dx d matrices, it is often possible to find an increasing 
sequence {nk)ke'N such that the sequence of the matrices A'f, = An^+i ■ ■ ■ An,,^^ satisfies 
them. This is what we do in a second paper, by considering a set of three 7x7 sparse 
(0, l)-matrices and all the sequences (A„)„gN with terms in this set, in order to apply 
the results of the present to the multifractal analysis of some Bernoulli convolution. 
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Introduction 

We associate to any sequence {An)n£^ of nonnegative dxd matrices, tlie sequence (Pn)nGN 
defined by 

(1) Pn := Al . . . An for n > 1, and Pq := I (identity matrix of order d). 

Under some liypotlieses on the matrices An, Theorem 11.11 gives some precisions about 

the orders of growth of the columns of P„ and the limits of their normalized columns 

PnV 

(we use the norm-sum). This theorem implies the convergence of ^^^"^^^ when ^ is a 
positive column vector (Corollary II. 2p . Our motivation is to apply this theorem - in [1] 



- in order to study some continuous singular measure defined by infinite convolution of 
discrete measures (see [T3] for the Bernoulli convolutions, and [TT] for their multifractal 
analysis by means of products of matrices). 

PnV 

Example 0.1. Let us give an example of positive matrices for which ^^^^^^ converges, 

Pn " 

and the normalized columns of Pn converge to the same limit, but diverges. We 



P 



n 
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first consider the sequence of nonnegative matrices (^'^^ g J associated to a sequence 

of positive integers (an)n.eN- We associate to (a„,)„gN the continued fraction (see [6] for 
instance): 



a 



ai, a2, ■ ■ ■ e N. 



ai 



0-2 



It is known that both normalized columns of the product matrix P, 



(a) 



converge to 



a 

a + 1 V 1 



ai 1 
1 



Consequently, for any nonnegative vector V, 



1 



(2) 



lim 

n— >oo 



a 

a + 1 V 1 



Let us prove that 



1/3(0!) I 



diverges excepted in case the an are eventually constant or tend 



to infinity: indeed if it converges, the ratio of the entries of its first row, let r„, either 
converges or tends to infinity; using the classical results about the continued fractions 



we have r„ 



Pn 

Pn-l 



> 1 and 



(3) 



1 



"n ' n 



rn-1 



so the integer an tends to a finite value (resp. to infinity) if Tn do 
not a product oj 



Notice that Pn°'^ is not a product of positive matrices, but Pg"'' is the product of the 



positive matrices 



0'2k 



for k 



n, and this gives the required 



counterexample: let us prove that 
(i) for any nonnegative vector V , 



converges to a limit that do not depend on V ; 



in particular the columns of Pj"^ converge to this limit; 



(U) 



2n 



P 



(a) 



2n 



diverges excepted if the a2„ are eventually constant or tend to infinity. 



Indeed (i) is an obvious consequence of and for proving (ii) we suppose that 



Die) 
2n 



P 



(a) 



2n 



converges. Then the ratio r2n of the entries of its first row tends to a finite value or to 
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1 



infinity. Now ^ implies = o-in 



1 



and 



0'2n-l 



(4) 



a2n < < a2n + C, whcrc C = sup 

k 



( \ 

1 



1 

r2k 



Consequently, 

_ in case ^ (finite) one has C < 1 and, for n large enough, \r\ < otherwise 

the fractional part of r2n should he close to 1 in contradiction with moreover, \r\ < 
f2n < [rj + 1 and a2n = ['"J for n large enough; 

_ in case r2n tends to infinity, a2n also do. 



The conditions for the convergence of the normahzed columns of P„ are also quite different 
from the conditions for the convergence of the matrix P„ itself (given for instance in 
[H El El in El El El [ID])- Let for instance Si, S2, ■ ■ ■ be a sequence of stochastic matrices; 
as noticed in [H Theorem 1] , if Si . . . Sn converges, the rows of the limit matrix are 
obviously left-eigenvectors - for the eigenvalue 1 - of each matrix S such that the set 
{n ; Sn = S} is infinite. Let (5'„)ngN be a sequence of stochastic 2x2 matrices, not 
diagonal and with positive diagonal entries, belonging to a finite set, such that Sn. has 

the form ^^_^^,xG]0, 1[ for infinitely many n, and Sn is not lower triangular for 

infinitely many n. The left-eigenvector of ^ ^ 1 x ^ eigenvalue 1 is ( 1 ) , 

which is not left-eigenvector of any non lower triangular matrix. So, Si . . . Sn diverges 
while - according to [17] - the normalized columns of 5*1 ... converge. 

In [12] and [17| we have treated the case of the products of 2 x 2 nonnegative matrices 
taken in a finite set; the conditions on this set, for the simple convergence, are quite 
different from the ones for the uniform convergence, that we use in [12] to study some 
Bernoulli convolutions. In [TT], [13] and [2] we also consider such measures, but only 
[Ti] requires to apply Theorem 11.11 because . in the case we consider, it is not possible to 
use matrices smaller than 7x7. Accordingly we think that Theorem 11.11 makes possible 
the study of the Bernoulli convolutions associated to the Pisot numbers, although the 
associated matrices may be large and far from satisfying the hypotheses of this theorem. 
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1. Properties of the columns of P„ under some hypotheses 



We use the following notation for the set of the indexes of the nonnuU entries in a column 

(Xi \ / mil • • • rnid 

'■ , or in a column of a matrix M = I : ' • . : 
Xd / \mdi ... rudd 

c{X) := {i ; Xij^ 0} and Cj{M) := {i ; ruij 7^ 0}. 

The first hypothesis we make means that the An = (^'^tj^^ echelon matrices in the 
following sense: 

(Hi): yn,j,j', Cj{An) C Cf{An) or Cj{An) ^ Cf{An). 

The second compare the entries of the j*^ and the j'^^ column of An. when Cj{An) contains 
strictly Cj'{An). 

{H2): Vn, z, J, j', if I e Cj (An) \ Cj, (An) then a^^ > 2 ^ . 

i' 

The third means that the ratio between two nonnuU entries located in a same column of 
a same matrix An is bounded independently on n: 



in) 



(H^) : \/n, i, i', j such that 4/ ^ 0, ^ < M < 00. 

The theorem below classify the columns of P„, with respect to their order of growth - this 
classification depends on n - and asserts the existence of an unique limit in direction for 
each order of growth. In the corollary we deduce the convergence of ^^^"^^^ when y is a 

positive vector. The norm we use is the norm-sum, and we call (£'j)i<j<d the canonical 
basis of the set of the d-dimensional column vectors. 

Theorem 1.1. For any sequence of nonnegative dxd matrices {An)neN that satisfy (Hi), 
{H2) and (i^^^), there exists some d-dimensional column vectors of norm 1, let Vi, . . . , Vs, 
and there exists for any n G N some nonempty and disjoint subsets of {l,...,d}, let 
D'('\...,Df\ such that 

3 e Dt^ = V. + o(l) and 3 ^ U^f ^ ^ P„E, = 0, 

1 1 n 1 1 

3 G =^ c (PnEj) = c{Vi) for n large enough, 
Ki' ^c{V,)Dc{Vi>), 

i < i' ^ lim I sup |^,^"^^''|,' | = 0. 

» jgD^"', i'sD'r' II^'^^jI 
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Corollary 1.2. Let be a sequence of nonnegative d x d matrices that satisfy 

(Hi), {H2) and (^H^^^ , and let V be a nonnegative vector such that PnV is not eventually 
null. Denoting by in{V) the smallest integer i such that c{V) fl D^"''' ^ % we have 



PnV ,r ^ ( 



Vji 

where max — means the greater ratio between the nonnull entries in V . In particular if 

PnV 

V is positive, lim exists and do not depend on V . 

Before proving Theorem II. we deduce directly from the hypotheses (Hi), {H2) and 
[H^] the following 

Lemma 1.3. // the matrices An satisfy (Hi), {H2) and (Z/^^), the products Pm,n = 
Am+i ...An= (^^17'"'') form<n satisfy 

{Hi): ym,n,j,j', Cj{P,n,n) ^ Cj>{Pm,n) OrCj{Pm,n) 5 Cj'{Pm,n); 

{H'2): ym,n,t,j,f, iftec,{Pm,n)\cAPm,n) ^/^en pj^'") > 2""'" || P„E,v || ; 

(m,n.) 

fH^Md\ . y ■ •/ ■ g^^i^ ^j^^^ , q ^ 2Mrf. 

\ / 1 1 1 ij jT'i J / ! (m,n) — 

Proof: Let us first prove {H2) by induction on n. It is true at the rank n = m + 1 because 
Pm,m+i = Am+1- Supposc it is true at some rank n — 1 > m + 1 and let be such 

that i e Cj{Pm,n) \ Cj'{Pm,n) that is, 

(5) ^ = 5:pSr"-'^4? and = pp^^ = Y.pt''"'^ af) . 



k k 



There exists such that p\^^ "'koj ^ ^^'^ Piko fcoj' ~ consequently a^^j 7^ and 
, = that is, ko G Cj{An)\cj'{An) and we can use the hypothesis (-^2): cik!^j > 2 a^^). 



(n) 



Moreover, from ([5]) any k such that a[."] 7^ satisfies p[^'^ = hence i G Cko{Pm,n-i) \ 
Ck{Pm,n-i) and - by the induction hypothesis 



(m,n-l) ^ 9n-l-m (m,n-l) 
Piko — ■ 2-^Pi'k 

Since p-^"'"^ > T>ik^~^^ ^^k^ij — ^ ^^pIa!^'""^^'^!:/ 5 where we sum for the indexes k such that 

(n) 
kj' 



k 

al") 7^ 0, we deduce 



(m,n) 
Pij 



k \ i' J 
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Pm.n satisfies the condition (Hi) because its columns are nonnegative linear combina- 
tion of the ones of A^+i and consequently each set Cj{Pm,n) is equal to some Cji{Am+i), 

fe{l,...,d}. 

It remains to prove that Pm,n satisfies (if^^^'^) . Let 



M^,.„. = max 



(m,n) 



(m,n) 

Pi'j 



that is, the maximum for E {1, . . . , d} such that p\7^'"'^ ^ 0. For such z, i', j we have 



(m,n-l) (n) 
(m,n) Pik '^kj 

Pij _ k 



{m,n) (m.,n— 1) (n) 

k '^kj 



{m,n) \- (r 
Pi'j 2^ Pi' 



k 



and, denoting by K the set of indexes k such that p-/*^'" ^ 0, 

E (m,n-l) (n) ST^ (m,n-l) (n) 
Pik (^kj l^Pik (^kj 

Pij keK , k<^K 



fe) = — + 

^ ' Jm.") (m,n-l) (n) ST^ (rn,n-\) in)- 

Pi' 3 2-^Pi'k ^kj 2-^Pi'k ^kj 

keK k£K 

In the first term we bound by In the second we use the inequality 

"^ptk'^'^^^tj - Pek'^~^^ ^tjj ' where ko is one of the indexes k such that 7^ 0. 

keK 

Now for k ^ K one has p\7l'^~^^ = while P-/™^"""'^^ 7^ 0, and {H2) implies that, for any 



^' plr'^"^'' — 2n-i-m P^^o^~^^ ' Using these bounds we deduce from ([6]) - after simplifica- 
tions - that 



(m,n) .i— ^ '^■J 



(m,n) - + 2"-l-™ H 

Pi'j ^koj 

< Mrr, n-1 + ; Md. 

The inequalities M^^n < ^m,n-i + ^— [z^^c^ forn = m+2,m+3,... , and Mm,m+i < 
prove that Mm,n < 2Md. ■ 

Remark 1.4. ('zj The condition (i/^^'^) m Lemma \1.3\ implies that the nonnull entries 
in Pm,n satisfy 

p-r' > 



(m,n) ^ ll -Pm,n-^j; I 



2Md2 • 

(a) Let j,j' be such that Cj{Pm^n) \ Cj'{Pm,n) contains at least one element, let iq. For any 
i G Cj{Pm,n) CLnd for any i' one has - from {H2) and (^H'^^^'^) 

i = i z£ <; (Ifld - ■ < 

(m,n) (m,n) (m,n) — On—m 11 p 7?. 1 1 — On—m 



ASYMPTOTIC PROPERTIES OF THE COLUMNS IN THE PRODUCTS OF NONNEGATIVE MATRICES 

Proof of the theorem: We can associate to any dxd nonnegative matrix M a matrix whose 

nonnull columns are the normahzed nonnuU columns of M: this is the matrix MD[M) 
/ Ai ... \ 



where D{M) :-- 



Aa ... 



1/ \\MEi\\ if \\MEi\\ ^ 
else. 



and Aj = 

V ... \d J 

This matrix belongs to the compact set 

/C := {M nonnegative dxd matrix ; Vz, HME'jll G {0, 1}} . 
Lemma 1.5. There exists an increasing sequence of integers, let {nk)keN, such that 

(i) the limit matrix P := lim Pn^D {Pn^) exists, and lim {uk+i — Uk) = oo; 

k—^oo k—^oo 

(ii) Cj{PnJ = Cj{P) for any j and k; 

\\P E -iW 

(Hi) lim — exists in the compact [0, oo] for any couple such that the j'*^ 

k-^co \\Pn„Ej\\ 

column of P differ from the null column. 

Proof: (i) In the compact set /C the sequence (P„-D(P„))„gN admits a convergent subse- 
quence, we take again a subsequence in order to have lim (n^+i — n^) = oo. 

k—^oo 

(ii) We replace {nk)km by some subsequence, in order that the set Cj{Pn^) is the same for 
any k. Hence - using (if^^*^) - it is also Cj{P). 

(iii) We replace again {nk)km by a suitable subsequence. ■ 

Let us denote by C the (may be empty) set of the indexes of the null columns of P, and 
by Ci, . . . , the equivalence classes defined from the relation 

J7^/^ hm ^^G]0,oo[, (j,/G{l,...,4\C), 

ordrered in such a way that lim "'^ = for j G Ca and j' G Ca+i- Let a{j) be the 
integer such that j G Ca{j). 

Lemma L6. (i) For any a G {1, ... ,7} the sets Cj{P) and Cj{Pn^^nk+i) do not depend on 
j G Cai we denote them by Ca{P) and Ca{Pnk,nk+i) respectively. 

(ii) The dxd matrix Pn^ (resp. Pnllnk+i) whose j^^ column is the one of Pn^. (resp. 
Pnk,nk+i) '^fj ^ Ca, and whose other columns are null, has a rectangular positive Ca{P) xCa 
(resp. Ca{Pnk,nk+i) ^ ^a) submatrix; its other entries are null. 

(iii) Denoting - for any a G {1, ... ,7} and k ^ N - by b = b{a, k) the smallest integer 
such that Ca{Pnk,nk+i) H Cfe 7^ 0, we have Cb{P) = Ca{P) and 

pi'^) ^ pW pi"-) 

fc — > 00 
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meaning that each nonnull entry of the left matrix is equivalent to the corresponding entry 
of the right one when A; — t- oo. 



Proof: (i) By contraposition let j, j' be such that Cj{P) is contains strictly Cj'{P). From 
Remark II .4( 11). ^| " || < has limit when n — )■ cxd hence a{j') > a{j). 

Suppose now that Cj{Pnf.^nk+i) contains strictly Cj/{Pn^.^nk+i)- Denoting - for any m and 
any n > m - by p^™'"'* the entries of Pm,n = ^m+i • • • ^n. one has 

(7) Pnk,.E, = J2l^^'''''^^P^^E, and Pn,,,E, = J^p^^^-^^ P^^^E,. 

i i 
^(nfc,nfc+i) 2Md 

According to Remark ll.4n i). — r r < • Denoting U = i 1 ... 1 ) one 

^ ^ ^' (nfc,nfc+i) — 2"fe+l-"fc V J 

Pij 



deduce 



\Pn,^^Ef\\ _ UPn,^,E^' ^ 2Md 



and a(j') > a(j). 

(ii) The entry of the i^^ row and j^^ column of P^^ (resp. Pnllnk+i) is positive if i G Ca{P) 
(resp. i G Ca(P„j.^„^.^J) and j G C^, otherwise it is null. 

(iii) We consider the first relation in ([7]) when j G Ca and we put b = b{a,k); we can 
restrict the summation to z G Ca{Pnk,nk+i) ■ It is clear that, for any i E Cb = b{a, k) and any 

"ilnmn vpntrir . j- n 

I J k 



i' G Cb' with b' > b, the entries of the column vector p^^*"^^^^^ Pnf,Eii are negligible when 



/c — )■ oo with regard to the ones of pl^'""'''^^^ Pn^Ei (we use the bound -A^ r < 2Md), 

Pij 

so the dominant terms in this summation are the ones with i G Ca{Pn^,nk+i) ^ ^h- 

Moreover, the same relation implies Ch{P) = Ca{P) because Cb'{P) C Cb{P) for any b' > b.u 

Let us prove that the column vectors PEj are equal for any j in a given class Ca- For this 
we can use the Birkhoff 's contraction coefficient defined in the third chapter of Seneta's 
book (see [I6] or [5]), after extending the required results of this chapter in the following 
way: 

Lemma 1.7. For any square or rectangular nonnegative matrix A = (aij) whose nonnull 
entries are the ones with indexes {i,j) in some product set, let I a x Ja, we put 



e{A)= max ^ and t{A) = ^^^. 
{i,j,k,i)eilxJl auajk V^(^) + 1 

For such matrices A and B, if AB exists and is not the null matrix we have 

\n{e{AB)) < ln{e (A)) -riB). 
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Proof: We can reduce the problem to the case where A and B are positive. Indeed, 
denoting by Aj j the submatrix constituted by the entries of A with indexes in some 
product set I x J, the hypothesis AB ^ implies 

Ja n Jb 7^ 0, Ic = Ia, Jc = Jb, Ci^ 
and, clearly, 9{AB) = 9{Cj^j^), ^(A.^.j^n/J < OiA) and r(5j^n/B,jJ < ^iB). 
From now A and B are positive; suppose first that i? is a square matrix. Then from 
[T6| Theorem 3.12], t{B) is the Birkfoff's contraction coefficient of B. Denoting by d the 
projective distance defined in [161 Section 3.1] we have 

MOiA)) = max d(Ri, R2) and \n(e(AB)) = max d(RiB, R2B) 

i?l,i?2 Rl,R2 

(maximums for any (i?i,i?2) couple of rows of A). Hence the inequality ln{6{AB)) < 
ln(^(y4)) ■ t{B) results from the following - which is an immediate consequence of the 
definition of the Birkfoff's contraction coefficient: 

d{RiB, R2B) < d{Ri, R2) ■ t{B). 

Suppose now i? is a m x n matrix with m < n. For any m x m submatrix B' of B we 
have 

\n{e{AB')) < \n{9{A)) -riB'). 

We deduce \n{6{AB)) < ln{6{A)) ■ t{B) because 0{B) and 9{AB) are the maximums - 
for B' m X m submatrix of i? - of 0{B') and 6[AB') respectively, and because r(-) is an 
increasing function of 6{-). 

Suppose B is amxn matrix with m > n and let B" he amxm supmatrix of B obtained 
by repeating some columns of B. Then \n{e{AB)) = \n{e{AB")) < \n{e{A)) ■ t{B") = 
He{A)) ■ t{B). m 

Lemma 1.8. The column vectors PEj for j in a same class Ca are equal. 

Proof: First, a straightforward property of 6 is that 0(C) < ( — ——] diD) if both 

\l-ej 

matrices C = {cij),D = (dij) are as in Lemma 11.71 with Ic x Jc = Id x Jd, and if 

1 — 5 < < 1 for any (i, j) ^ Ic x Jc- According to Lemma [LOT iii) . these conditions 
dij 

are satisfied by C = Pi^+i and D = Pn^ Pnk,nk+^ for k large enough hence - using Lemma 
Owith A = pS and B = p/"^ 



Now for any a G {1, . . . ,7}, ^ {Pn^^ tends to 9 (p(")) when k — )■ 00, where P^"-* is the 
dxd matrix which has the same j^^ column as P for j ^ Ca, and whose other columns are 
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null. From {H^^^^), 6 (PrSk+i) is bounded by {2Md)^ hence r (pi"k+i) is bounded by 
a constant p < 1, and ([8]) implies 

max In (0 (p("))) < 21n|^^^^ +p-maxln(^ (p("))) 

for any e > that is, max In (6* (P^"^)) = and each P*^"^ has rank 1. ■ 

It remains to extend these results concerning the columns of the P„j., to the ones of P„ 
for any n. The number of classes Ci, . . . - defined from the orders of growth of the 
columns of P„j. - may obviously depend on the choice of the sequence {nk)keN (see Remark 
11.101 below) so we must define some new classes by grouping together the C^. Clearly, 
there exists = < ai < ■ ■ ■ < = 7 such that the Vj := PEj are the same for 
j eDi:= [j Ca, and Vj ^ Vf if j G A and f G A+i- 

Lemma 1.9. Given z G {1, . . . , 6}, for any a G]aj_i, Oj] and for k large enough the integer 
b{a,k) defined in Lemma \1 . 0( Hi) belongs to ]aj_i,ai]. 

Proof: Let us first prove that Vf,(^a,k) = Va for k large enough. Let i G Cq, we have from 
Lemma ll.6( iii) 

p(a) 771 pW p{«) 771 



p(«) r;. 



fc — > 00 



pCf-) p(«) p 



where the l.h.s. is a normalized column of Pi^+i and the r.h.s. is a linear combination of 
the normalized columns of Pn^', the sum of the coefficients of this linear combination is 1. 
Making k 00 along a sequence of integers k such that b{a, k) has a given value, let b, 
we deduce Va = H- In other words, Va = Vb{a,k) excepted for a finite number of values 
of A;. 

Now for a G]aj_i,aj] and a' G]aj,ai+i], b{a,k) and b{a',k) cannot be in a same interval 
]aj,aj^i] for infinitely many k: if they do, Vi,(^a,k) = H(a',fc) hence Va = Va'. They are in 
two intervals, let ] 1 and 1 ] respectively. We cannot have j' < j: if we have. 

Lemma ll.6( in) implies that the order of growth of Prl^+i is less than the one of Pt'l , m 
contradiction with a < a'. 

Finally, let us choose one element a in each interval ]aj_i,aj]; for k large enough, the 
corresponding element b{a, k) belong to some interval jctji-i, ctjj with 1 < Ji < J2 < ■ ■ ■ < 
is < 5- This implies obviously ji = i for any i. ■ 

We can now define for any n > rii the set -Dj"^ involved in the theorem: let k = k{ri) 
be the integer such that < n < Uk+i and -D^-"^ the set of the indexes j such that 
Cj{Pn^^n) intersects Di and do not intersect Dii for i' < i. Let us first prove that D^^^ 
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is non empty for n large enough. Let j G -Dj and k = k{n); using Lemma II. 9[ for k 
large enough Cj{Pn,.,nk+i) intersects Di and do not intersect Di/ for i' < i. Since the j^^ 
column of Pnk,nk+i is a linear combination of the ones of Pnk,n , there exists j' such that 
Cj{Pnk,nk+i) = <^j'{fnk,n) and Consequently j' G -D,/"''. It remains modify the set D^"^ for 
a finite number of integers n, in order that D^"^ ^ holds for any n (the assertions of 
Theorem 11.11 do not depend on this modification). 

The column of P„ is PnEj = Pn^Puk^nEj] it is a linear combination of the columns 
of Pn,,, and the ratio between two coefficients of this linear combination is bounded from 
Lemma 11.31 Using the different orders of growth of the columns of P^^. , we deduce that 
the first assertion in Theorem 11.11 is true - provides we call now Vi the common value of 
Va for a G]aj_i,aj]. The second results from Lemma [LST ii). and the two last assertions 
from the definition of the classes Ca- ■ 

Proof of Corollary \ 1.2: We denote by i the integer iniV) involved in Corollary 11.21 we 
have c{V) n D^"^ ^ and c{V) n D^^^ = for z' < z. From Theorem [TH 



PnV = VjPnE,+ J2 ^fPnEj 

(9) 



Y \\PnEj\\ V, + o (\\PnV\\ ■ max^ j . 



Making the product by ( 1 ... 1 ) one deduce 

(10) ||P„l^|| = Y \\PnEj\\ + o (\\PnV\\ ■ max 



and, from fl9|) and (JTOj) . 



P^V = \\PnV\\ Vi + o{ \\PnV\\ ■ max^ 



PnV 

which proves the corollary. In particular if V is positive, Vi = Vi is the limit of 

1 1 -Tfi, V I 



Example 1.10. Let us see on a trivial example that the classes D^"^^ of Theorem \1.1\ - 
related to the different limit vectors - are not the classes Ca - related to the different orders 
of growth in the choosen subsequence (P^^J^eN- Let {nk)km be an increasing sequence of 
integers such that hk = Uk+i — Uk tends to oo when k — )■ oo, and {An)n&n ^ sequence of 
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matrices defined by 



1 1 


2 


1 


1 \ 




( 36 








\ 


2 


1 


1 


1 


1 — 


30 


6 




















30 


4 


2 













^) 




I 30 


4 


1 










/ 36 


















An 


36 


6 
















36 
36 


6 
6 



1 




1/ 



/or n ^ {1, ni + 1, ^2 + 1, . . . }, 



for n e {ni + 1, n2 + 1, . . . }. 



Putting Sfc = ^ — andak = ^ 

/ (l + 2Sfc + Sfc2 

(2 + Sfc + 



j=0 



1 

36" 



(^wzi/i 77-0 = 0^ one has 



V 



^fe)36" 
afc)36"^-i 





- (2sfc + 2)6" 
(2sfc + 1)6" 



2)6"fc-i - 2^" 2'"= - 1 

ijg"*-! - 2^* 2^*^ - 1 





1 \ 

1 


0/ 



hence there exists - for this subsequence {Pnjkm - four orders of growth and four classes Ca- 
Nevertheless 

/ (1 + 2sk + + afe)36"^ {2sk + 2)6"^' 
(2 + Sfc + + afc)36"^ (2sfc + 1)6"'= 





v 



1 1\ 

1 1 


0/ 

and there exists three orders of growth for {Pnk+i)km- The classes involved in Theorem \l.l\ 
are = {1},D^^ = {2} and D^'' = {3,4}. The normalized columns of Pn have three 
distinct limits. 

Question 1.11. For instance what are the necessary and sufficient conditions on the 
triangular by blocks matrices An = ( W ' ' ^^^^^ ^"' ^" ^" positive, for 



Ai...AnV 
\Ai...AnV\\ 



to converge for any positive column vector V ? 
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